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Abstract 

Evolutionary dynamics have been traditionally studied in the context of homogeneous populations, 
mainly described by the Moran process Recently, this approach has been generalized in [TB] by 

arranging individuals on the nodes of a network (in general, directed). In this setting, the existence of 
directed arcs enables the simulation of extreme phenomena, where the fixation probability of a randomly 
placed mutant (i.e. the probability that the offsprings of the mutant eventually spread over the whole 
population) is arbitrarily small or large. On the other hand, undirected networks (i.e. undirected graphs) 
seem to have a smoother behavior, and thus it is more challenging to find suppressors/amplifiers of 
selection, that is, graphs with smaller/greater fixation probability than the complete graph (i.e. the 
homogeneous population). In this paper we focus on undirected graphs. We present the first class of 
undirected graphs which act as suppressors of selection, by achieving a fixation probability that is at most 
one half of that of the complete graph, as the number of vertices increases. Moreover, we provide some 
generic upper and lower bounds for the fixation probability of general undirected graphs. As our main 
contribution, we introduce the natural alternative of the model proposed in [T6'. In our new evolutionary 
model, all individuals interact simultaneously and the result is a compromise between aggressive and 
non-aggressive individuals. That is, the behavior of the individuals in our new model and in the model 
of [16] can be interpreted as an "aggregation" vs. an "all-or-nothing" strategy, respectively. We prove 
that our new model of mutual influences admits a potential function, which guarantees the convergence 
of the system for any graph topology and any initial fitness vector of the individuals. Furthermore, 
we prove fast convergence to the stable state for the case of the complete graph, as well as we provide 
almost tight bounds on the limit fitness of the individuals. Apart from being important on its own, this 
new evolutionary model appears to be useful also in the abstract modeling of control mechanisms over 
invading populations in networks. We demonstrate this by introducing and analyzing two alternative 
control approaches, for which we bound the time needed to stabilize to the "healthy" state of the system. 

Keywords: Evolutionary dynamics, undirected graphs, fixation probability, potential function, 
Markov chain, fitness, population structure. 



1 Introduction 

Evolutionary dynamics have been well studied (see Ell71l5]|14ll22|j24ll25j ). mainly in the context of homogeneous 
populations, described by the Moran process ,18.,20) . In addition, population dynamics have been extensively 
studied also from the perspective of the strategic interaction in evolutionary game theory, cf. for instance |10f - 
UnillS]. One of the main targets of evolutionary game theory is evolutionary dynamics (see [UlIlS]). Such 
dynamics usually examine the propagation of intruders with a given fitness to a population, whose initial 
members (resident individuals) have a different fitness. In fact, "evolutionary stability" is the case where no 
dissident behaviour can invade and dominate the population. The evolutionary models and the dynamics 
we consider here belong to this framework. In addition, however, we consider structured populations (i.e. in 
the form of an undirected graph) and we study how the underlying graph structure affects the evolutionary 
dynamics. We study in this paper two kinds of evolutionary dynamics. Namely, the "all or nothing" case 
(where either the intruder overtakes the whole graph or die out) and the "aggregation" case (more similar 
in spirit to classical evolutionary game theory, where the intruder's fitness aggregates with the population 
fitness and generates eventually a homogeneous crowd with a new fitness). 

In a recent article, Lieberman, Haucrt, and Nowak proposed a generalization of the Moran process by ar- 
ranging individuals on a connected network (i.e. graph) [16j (see also [21j). In this model, vertices correspond 
to individuals of the population and weighted edges represent the reproductive rates between the adjacent 
vertices. That is, the population structure is translated into a network (i.e. graph) structure. Furthermore, 
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individuals (i.e. vertices) are partitioned into two types: aggressive and non-aggressive. The degree of (rel- 
ative) aggressiveness of an individual is measured by its relative fitness; in particular, non-aggressive and 
aggressive individuals are assumed to have relative fitness 1 and r > 1, respectively. This modeling approach 
initiates an ambitious direction of interdisciplinary research, which combines classical aspects of computer 
science (such as combinatorial structures and complex network topologies), probabilistic calculus (discrete 
Markov chains), and fundamental aspects of evolutionary game theory (such as evolutionary dynamics). 

In the model of iTB], one mutant (or invader) with relative fitness r > 1 is introduced into a given 
population of resident individuals, each of whom having relative fitness 1. For simplicity, a vertex of the 
graph that is occupied by a mutant will be referred to as black, while the rest of the vertices will be referred 
to as white. At each time step, an individual is chosen for reproduction with a probability proportional to 
its fitness, while its offspring replaces a randomly chosen neighboring individual in the population. Once u 
has been selected for reproduction, the probability that vertex u places its offspring into position v is given 
by the weight Wuv of the directed arc (uv) . This process stops when either all vertices of the graph become 
black (resulting to a fixation of the graph) or they all become white (resulting to extinction of the mutants) . 
Several similar models have been previously studied, describing for instance influence propagation in social 
networks (such as the decreasing cascade model [miS]), dynamic monopolies [3], particle interactions (such 
as the voter model, the antivoter model, and the exclusion process [llElITT]), etc. However, the dynamics 
emerging from these models do not consider different fitnesses for the individuals. 

The fixation probability fa of a graph G — (V, E) is the probability that eventually fixation occurs, i.e. the 
probability that an initially introduced mutant, placed uniformly at random on a vertex of G, eventually 
spreads over the whole population V, replacing all resident individuals. One of the main characteristics 
in this model is that at every iteration of the process, a "battle" takes place between aggressive and non- 
aggressive individuals, while the process stabilizes only when one of the two teams of individuals takes over the 
whole population. This kind of behavior of the individuals can be interpreted as an all-or-nothing strategy. 

Lieberman et al. [H] proved that the fixation probability for every symmetric directed graph (i.e. when 
Wuv = Wyu for every u, v) is equal to that of the complete graph (i.e. the homogeneous population of the Moran 
process), which tends to 1 — ^ as the size n of the population grows. Moreover, exploiting vertices with zero 
in-degree or zero out-degree ("upstream" and "downstream" populations, respectively), they provided several 
examples of directed graphs with arbitrarily small and arbitrarily large fixation probability |16j . Furthermore, 
the existence of directions on the arcs leads to examples where neither fixation nor extinction is possible (e.g. a 
graph with two sources). 

In contrast, general undirected graphs (i.e. when {uv) £ E li and only if {vu) € E for every u, v) appear to 
have a smoother behavior, as the above process eventually reaches fixation or extinction with probability 1. 
Furthermore, the coexistence of both directions at every edge in an undirected graph seems to make it more 
difficult to find suppressors or amplifiers of selection (i.e. graphs with smaller or greater fixation probability 
than the complete graph, respectively), or even to derive non-trivial upper and lower bound for the fixation 
probability on general undirected graphs. This is the main reason why only little progress has been done so 
far in this direction and why most of the recent work focuses mainly on the exact or numerical computation 
of the fixation probability for very special cases of undirected graphs, e.g. the star and the path [^{B]. 

Our contribution. In this paper we overcome this difficulty for undirected graphs and we provide the first 
class of undirected graphs that act as suppressors of selection in the model of [16] , as the number of vertices 
increases. This is a very simple class of graphs (called clique-wheels) , where each member Gn has a clique 
of size n > 3 and an induced cycle of the same size n with a perfect matching between them. We prove 
that, when the mutant is introduced to a clique vertex of G„, then the probability of fixation tends to zero 
as n grows. Furthermore, we prove that, when the mutant is introduced to a cycle vertex of G„, then the 
probability of fixation is at most 1 — ^ as n grows (i.e. to the same value with the homogeneous population 
of the Moran process). Therefore, since the clique and the cycle have the same number n of vertices in G„, 
the fixation probability /g„ of G„ is at most ^(1 — i) as n increases, i.e. G„ is a suppressor of selection. 
Furthermore, we provide for the model of 1161 the first non-trivial upper and lower bounds for the fixation 
probability in general undirected graphs. In particular, we first provide a generic upper bound depending on 
the degrees of some local neighborhood. Second, we present another upper and lower bound, depending on 
the ratio between the minimum and the maximum degree of the vertices. 

As our main contribution, we introduce in this paper the natural alternative of the all-or-nothing ap- 
proach of [12], which can be interpreted as an aggregation strategy. In this aggregation model, all individuals 
interact simultaneously and the result is a compromise between the aggressive and non-aggressive individu- 
als. Both these two alternative models for evolutionary dynamics coexist in several domains of interaction 
between individuals, e.g. in society (dictatorship vs. democracy, war vs. negotiation) and biology (natural se- 
lection vs. mutation of species). In particular, another motivation for our models comes from biological 
networks, in which the interacting individuals (vertices) correspond to cells of an organ and advantageous 
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mutants correspond to viral cells or cancer. Regarding the proposed model of mutual influences, we first 
prove that it admits a potential function. This potential function guarantees that for any graph topology 
and any initial fitness vector, the system converges to a stable state, where all individuals have the same 
fitness. Furthermore, we analyze the telescopic behavior of this model for the complete graph. In particular, 
we prove fast convergence to the stable state, as well as we provide almost tight bounds on the limit fitness 
of the individuals. 

Apart from being important on its own, this new evolutionary model enables also the abstract modeling 
of new control mechanisms over invading populations in networks. We demonstrate this by introducing and 
analyzing the behavior of two alternative control approaches. In both scenarios we periodically modify the 
fitness of a small fraction of individuals in the current population, which is arranged on a complete graph 
with n vertices. In the first scenario, we proceed in phases. Namely, after each modification, we let the 
system stabilize before we perform the next modification. In the second scenario, we modify the fitness of a 
small fraction of individuals at each step. In both alternatives, we stop performing these modifications of the 
population whenever the fitness of every individual becomes sufficiently close to 1 (which is considered to be 
the "healthy" state of the system). For the first scenario, we prove that the number of phases needed for the 
system to stabilize in the healthy state is logarithmic in r — 1 and independent of n. For the second scenario, 
we prove that the number of iterations needed for the system to stabilize in the healthy state is linear in n 
and proportional to rln(r — 1). 

Notation. In an undirected graph G — {V,E), the edge between vertices u G V and v G V is denoted 
by uv £ E, and in this case u and v are said to be adjacent in G. If the graph G is directed, we denote 
by (uv) the arc from u to v. For every vertex m e in an undirected graph G = (V^E), we denote by 
N{u) = {v gV \ uv £ E} the set of neighbors of u in G and by deg(u) — |iV(w)|. Furthermore, for any A; > 1, 
we denote for simplicity [fc] = {1, 2, . . . , fc}. 

Organization of the paper. We discuss in Section [2] the two alternative models for evolutionary dynamics 
on graphs. In particular, we formally present in Section 12.11 the model of |16| and then we introduce in 
Section f2.2l our new model of mutual influences. In Section[3]we first provide generic upper and lower bounds 
of the fixation probability in the model of [16j for arbitrary undirected graphs. Then we present in Section f3.3l 
the first class of undirected graphs which act as suppressors of selection in the model of [TS] , as the number of 
vertices increases. In Section|4]we analyze our new evolutionary model of mutual infiuences. In particular, we 
first prove in Section 14.11 the convergence of the model by using a potential function, and then we analyze in 
Section the telescopic behavior of this model for the case of a complete graph. In Section[S]we demonstrate 
the use of our new model in analyzing the behavior of two alternative invasion control mechanisms. Finally, 
we discuss the presented results and further research in Section El 

2 All-or-nothing vs. aggregation 

In this section we formally define the model of [TB] for undirected graphs and we introduce our new model of 
mutual influences. Similarly to jl6j . we assume for every edge uv of an undirected graph that Wuv = ^cgu ^^'^ 
Wyu = , i-e. once a vertex u has been chosen for reproduction, it chooses one of its neighbors uniformly 
at random. 

2.1 The model of Lieberman, Hauert, and Nowak (an all-or-nothing approach) 

Let G = (y, E) be a connected undirected graph with n vertices. Then, the stochastic process defined in [16] 
can be described by a Markov chain with state space S = 2^ (i.e. the set of all subsets of V) and transition 
probability matrix P, where for any two states Si, S2 C V, 

if 82 = SiU {v} and v ^ Si 
if 5i = 5*2 U {v} and v S2 

(1) 

if S2 — Si 
otherwise 

Notice that in the above Markov chain there are two absorbing states, namely and V, which describe 
the cases where the vertices of G are all white or all black, respectively. Since G is connected, the above 
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Markov chain will eventually reach one of these two absorbing states with probability 1. If we denote by hy 
the probability of absorption at state given that we start with a single mutant placed initially on vertex 
w, then by definition fa — . Generalizing this notation, let hs be the probability of absorption at V 

given that we start at state S ^V, and let h = [hs]s<zv- Then, it follows that vector h is the unique solution 
of the linear system h = P ■ h with boundary conditions /ig = and hy = 1 ■ 

However, observe that the state space 5 = 2^ of this Markov chain has size 2", i.e. the matrix P — [Psi^S^] 
in H]) has dimension 2" x 2". This indicates that the problem of computing the fixation probability fa of a 
given graph G is hard, as also mentioned in 16^. This is the main reason why, to the best of our knowledge, all 
known results so far regarding the computation of the fixation probability of undirected graphs are restricted 
to regular graphs, stars, and paths [i Ml]16ll21) . In particular, for the case of regular graphs, the above Markov 
chain is equivalent to a birth-death process with n — 1 transient (non-absorbing) states, where the forward 
bias at every state (i.e. the ratio of the forward probability over the backward probability) is equal to r. In 
this case, the fixation probability is equal to 

1 1 _ 1 

(2) 



cf. |21| . chapter 8. It is worth mentioning that, even for the case of paths, there is no known exact or 
approximate formula for the fixation probability [S]. 

2.2 An evolutionary model of mutual influences (an aggregation approach) 

The evolutionary model of [T^ constitutes a sequential process, in every step of which only two individuals 
interact and the process eventually reaches one of two extreme states. However, in many evolutionary 
processes, all individuals may interact simultaneously at each time step, while some individuals have greater 
influence to the rest of the population than others. This observation leads naturally to the following model 
for evolution on graphs, which can be thought as a smooth version of the model presented in |16| . 

Consider a population of size n and a portion a G [0, 1] of newly introduced mutants with relative fitness r. 
The topology of the population is given in general by a directed graph G = {V,E) with \V\ = n vertices, 
where the directed arcs of E describe the allowed interactions between the individuals. At each time step, 
every individual u &V oi the population influences every individual w G V^, for which {uv) € E, while the 
degree of this influence is proportional to the fitness of u and to the weight Wuv of the arc {uv) . Note that 
we can assume without loss of generality that the weights Wuv on the arcs are normalized, i.e. for every fixed 
vertex u & V \t holds X](Mt))e£; "^^v — 1 ■ Although this model can be defined in general for directed graphs 
with arbitrary arc weights Wuvi we will focus in the following to the case where G is an undirected graph 
(i.e. {uiUj) G E' if and only if {ujUi) G E, for every and Wuv = dcg{u) ^'-'^ edges uv e E. 

Formally, let V ~ {ui,U2, ■ ■ ■ , Un} be the set of vertices and r„. (fc) be the fitness of the vertex Ui G V a.t 
iteration k > 0. Let S(fc) denote the sum of the fitnesses of all vertices at iteration k, i.e. S(fc) = X]r=i ^ui{k)- 
Then the vector r(fc -|- 1) with the fitnesses (fc -I- 1) of the vertices Ui G V at the next iteration fc -|- 1 is 
given by 

[ru, (fc -f 1), (fc + 1), . . . , r„„ (fc + 1)]^ = P ■ [r„, (fc), (fc), . . . , r„„ (fc)]^ (3) 

i.e. 

r(fc + 1) = F • r(fc) (4) 

In the latter equation, the elements of the square matrix P — [Pij]i j^i depend on the iteration fc and they 
are given as follows: 

{r (k) 
dcg(n,)S(fc) ' if ' ^ J' and u,uj e E 

0, if « 7^ j and UiUj ^ E (5) 

Note by (jl]) and (O that after the first iteration, the fitness of every individual in our new evolutionary 
model of mutual influences equals the expected fitness of this individual in the model of [IH] (cf. Section [O]) . 
However, this correlation of the two models is not maintained in the next iterations and the two models 
behave differently as the processes evolve. 

In particular, in the case where G is the complete graph, i.e. deg(ui) = n — I for every vertex Ui, the 
matrix P becomes 
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(6) 



(n-l)S(fc) (n-l)E(fe) " " (n-l)S(fc) 

The system given by (HI and ([5]) can be defined for every initial fitness vector r(0). However, in the case 
where there is initially a portion a G [0, 1] of newly introduced mutants with relative fitness r, the initial 
condition r(0) of the system in is a vector with an entries equal to r and with (1 — a)n entries equal to 1. 

Observation 1 Note that the recursive equation ^ is a non-linear equation on the fitness values r^^ (k) of 
the vertices at iteration k. 

Since by ([5]) the sum of every row of the matrix P equals to one, the fitness . (k) of vertex Ui after the 
(fc + l)-th iteration of the process is a convex combination of the fitnesses of the neighbors of Ui after the 
fc-th iteration. Therefore, in particular, the fitness of every vertex Ui at every iteration fc > lies between 
the smallest and the greatest initial fitness of the vertices, as the next observation states. 

Observation 2 Let rmin and rmax be the smallest and the greatest initial fitness in r(0), respectively. Then 
rmin < ^ui(k) < Tniax for every Ui £V and every k >0. 



Degree of influence. Suppose that initially an mutants (for some a G [0,1]) with relative fitness r > 1 
are introduced in graph G on a subset S C V oi its vertices. Then, as we prove in Theorem [5l after a certain 



number of iterations the fitness vector r{k) converges to a vector 



' I ' I • ■ 



'■g]^, for some value Tq . This 



limit fitness depends in general on the initial relative fitness r of the mutants, on their initial number an, 
as well as on their initial position on the vertices oi S C V. The relative fitness r of the initially introduced 
mutants can be thought as having the "black" color, while the initial fitness of all the other vertices can be 
thought as having the "white" color. Then, the limit fitness Tq can be thought as the "degree of gray color" 
that all the vertices obtain after sufficiently many iterations, given that the mutants are initially placed at the 
vertices of S. In the case where the an mutants are initially placed with uniform probability to the vertices 
of G, we can define the limit fitness tq of G as 



SCV, |S|=Qn 
van/ 



(7) 



For a given initial value of r, the bigger is tq the stronger is the effect of natural selection in G. 

Since Tq is a convex combination of r and 1, there exists a value fcsi^) € such that 

'^0 = /g,s(^) ■ r + {1 — fG,s{r)) ■ 1- Then, the value fo-sif) is the degree of influence of the graph G, given 
that the mutants are initially placed at the vertices of S*. In the case where the mutants are initially placed 
with uniform probability at the vertices of G, we can define the degree of influence of G as 

E fcsir) 

, S<ZV, \S\=an . 

lG{r) = — (8) 



Van, 



) 



Number of iterations to stability. For some graphs G, the fitness vector r{k) reaches exactly the limit 
fitness vector [ro,ro, . . . ,ro]"^ (for instance, the complete graph with two vertices and one mutant not only 
reaches this limit in exactly one iteration, but also the degree of influence is exactly the flxation probability 
of this simple graph). However, for other graphs G the fltness vector r{k) converges to [ro, ro, . . . , tq]^ 
(cf. Theorem [S] below), but it never becomes equal to it. In the first case, one can compute (exactly or 
approximately) the number of iterations needed to reach the limit fitness vector. In the second case, given an 
arbitrary e > 0, one can compute the number of iterations needed to come e-close to the limit fitness vector. 



3 Analysis of the all-or-nothing model 

In this section we present analytic results on the evolutionary model of |16j , which is based on the sequential 
interaction among the individuals. In particular, we first present non-trivial upper and lower bounds for the 
fixation probability, depending on the degrees of vertices. Then we present the first class of undirected graphs 
that act as suppressors of selection in the model of [16 , as the number of vertices increases. 
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Recall by the preamble of Section [2.2l that. similarly to [16j . we assumed that Wuv — ^^gu ^'^^ ^'"f ~ dcgv 
for every edge uv of an undirected graph G = {V,E). It is easy to see that this formulation is equivalent 
to assigning to every edge e = uv € E the weight We = Wuv = Wvu = 1, since also in this case, once a 
vertex u has been chosen for reproduction, it chooses one of its neighbors uniformly at random. A natural 
generalization of this weight assignment is to consider G as a complete graph, where every edge e in the clique 
is assigned a non- negative weight > 0, and We is not necessarily an integer. Note that, whenever = 0, 
it is as if the edge e is not present in G. Then, once a vertex u has been chosen for reproduction, u chooses 
any other vertex v with probability ^ "'"''^ — . 

Note that, if we do not impose any additional constraint on the weights, we can simulate multigraphs by 
just setting the weight of an edge to be equal to the multiplicity of this edge. Furthermore, we can construct 
graphs with arbitrary small fixation probability. For instance, consider an undirected star with n leaves, 
where one of the edges has weight an arbitrary small e > and all the other edges have weight 1. Then, the 
leaf that is incident to the edge with weight e acts as a source in the graph as £ -> 0. Thus, the only chance 
to reach fixation is when we initially place the mutant at the source, i.e. the fixation probability of this graph 
tends to as e — 0. Therefore, it seems that the difficulty to construct strong suppressors lies in the fact 
that unweighted undirected graphs can not simulate sources. For this reason, we consider in the remainder 
of this paper only unweighted undirected graphs. 

3.1 A generic upper bound approach 

In the next theorem we provide a generic upper bound of the fixation probability of undirected graphs, 
depending on the degrees of the vertices in some local neighborhood. 

Theorem 1 Let G = {V, E) he an undirected graph. For any uv G E, let Qu — 'YlixeN(u) dcgz '^'^'^ ~ 

^xeN(u)\{v} + HxeN(v)\{ii} T'^'i'^ 

fG < I . q^q Z^ I (9) 



rQu 



Proof. For the proof we construct a simple Markov chain Ai , in which the probability of reaching one of its 
absorbing states is at least the probability of fixation in the original Markov chain. Then, in order to provide 
an upper bound of the fixation probability in the original Markov chain, we provide an upper bound on the 
probability of reaching one of the absorbing states in Ai . 

Let w be a vertex that maximizes the probability of fixation, that is, u € argmax„gy Furthermore, 
assume that we end the process in favor of the black vertices when the corresponding Markov chain describing 
the model of |16) reaches three black vertices. To favor fixation even more, since u maximizes h^, we assume 
that, whenever we reach two black vertices and a backward step happens (i.e., a step that reduces the number 
of black vertices), then we backtrack to state u (even if vertex u was the one that became white). Finally, 
given that we start at vertex u and we increase the number of black vertices by one, we assume that we make 
black the neighbor w of u that maximizes the forward bias of the state {u, v}. Imposing these constraints (and 
eliminating self loops), we get a Markov chain A4, shown in Figure [T] that dominates the original Markov 
chain. That is, the probability that Ai reaches the state of three black vertices, given that we start at u, is 
an upper bound of the fixation probability of G. 

Pl , P2 
^ U, V J 

92 ^— ^ 

Figure 1: The Markov chain Ai. 
For the Markov chain Ai, we have that 

'Ylix&Niu) desx def Qu -, 
' + Z^a;GJV(u) dega: ^ ^" 

where N{u) is the set of neighbors of u. Also, 

ExeJV(u)\{«} dSi X {v)\{u\ dcg a: 
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Let now /i„ (resp. h^v) denote the probability of reaching 3 blacks, starting from u (resp. starting from the 
state {u,v}) in A4. We have that 

K = PiKv = Pi{p2 + q2hu) 

^ _ P1P2 _ ^ 

1 - piq2 r2 + rQu + ^ S"'^"" 1 ~ + rQu + 

dog u dcg u 

This completes the proof of the theorem. ■ 

Consider for instance a bipartite graph G = {U,V,E), where degw = di for every vertex u £ U and 
degi; = d2 for every vertex v G V. Then any edge of E has one vertex in U and one vertex in V. Using 
the above notation, consider now an arbitrary edge uv £ E, where u £ U and v G V. Then Qu = ^ and 
Quv = ^d"'" + ■ The right side of ([9]) is maximized when di < d2, and thus in this case Theorem [1] implies 
that fa < Ti — /di-i — d2-i\ ■ particular, for the star graph with n + 1 vertices, we have di — 1 and 

d2 — n. But, as shown in [16[ . the fixation probability of the star is asymptotically equal to 1 — ^, whereas 

2 1 

the above bound gives fstar < ^^^^ ^ ^ „_i = 1 — ^.^^-^Yo(i) ; which is quite tight. 
3.2 Upper and lower bounds depending on degrees 

In the following theorem we provide upper and lower bounds of the fixation probability of undirected graphs, 
depending on the minimum and the maximum degree of the vertices. 

Theorem 2 Let G = {V, E) he an undirected graph, where < A whenever uv G E. Then, the fixation 

probability Jq 0/ G, when the fitness of the mutant is r, is upper (resp. lower) bounded by the fixation 
probability of the clique for mutant fitness ri = rX (resp. for mutant fitness r2 = j). That is, 

Proof. For an arbitrary state S C V oi the Markov Chain (that corresponds to the set of black vertices 
in that state), let p+{S) (resp. p-{S)) denote the probability that the number of black vertices increases 
(resp. decreases). In the case where G is the clique, the forward bias at state S is equal to y^^^ = r, for 
every state S [TnHH]. Let Cs = {uv € E \ u € S, v ^ S} he the set of edges with one vertex in S and one 
vertex in V \ S. Then, 

P+^^^^^ ^n-\S\+r\S\degiu) ^^^^ 

{uveE I ues, v(^s} ' ' ' ' ' 

and 

P-^^^^ E iJl^ KM \ ^ (12) 

iu.eE lues. ^^^^n-\S\+r\S\deg{v) 

Now, since by assumption < A whenever uv G E, it follows that 

A E deg(w) ~ E deg(w) ~ ^ ^ deg(w) ^^'^'^ 

{iiveE I ues. vfs} ^ ' {uveE \ ues, vfs} ' {uveE \ ues, v<^s} ' 

By HH), HH), and ^ we get the following upper and lower bounds for the forward bias at state S. 

{<^,<rX (14) 
A p_(5) 

Notice that the upper and lower bounds of ([T^ for the forward bias at state S are independent of S. Therefore, 
the process stochastically dominates a birth-death process with forward bias j, while it is stochastically 
dominated by a birth-death process with forward bias rX (cf. equation ([2])). This completes the proof of the 
theorem. ■ 
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3.3 The undirected suppressor 

In this section we provide the first class of undirected graphs (which we call clique-wheels) that act as sup- 
pressors of selection as the number of vertices increases. In particular, we prove that the fixation probability 
of the members of this class is at most 5(1—7), i-e. the half of the fixation probability of the complete graph, 
as n — >■ 00. An example of a clique-wheel graph G„ is depicted in Figure 2(a) This graph consists of a clique 
of size n > 3 and an induced cycle of the same size n with a perfect matching between them. We will refer 
in the following to the vertices of the inner clique as clique vertices and to the vertices of the outer cycle as 
ring vertices. 



G„ 




71 A 



72 il 




(a) 



F3. 



73 ik 



(b) 



^-y^ — 1 

^ log ' n 



an — 1 

log ' n 



Sn_, 



Figure 2: (a) The clique- wheel graph Gn and (b) the state graph of a relaxed Markov chain for computing 
an upper bound of hi = hcUque- 



Denote by hcUque (resp. hring) the probability that all the vertices of G„ become black, given that we 
start with one black clique vertex (resp. with one black ring vertex). We first provide in the next lemma an 
upper bound on hcUque- 



Lemma 1 For any r E (l, |), 



hcii. 



< 



6n I 



3r 



1) 



O I - 

n 



Proof. Denote by Sk the state, in which exactly fc > clique vertices are black and all ring vertices are 
white. Note that 5*0 is the empty state. Denote by Fk the state where at least one ring vertex of On and 
exactly fc > clique vertices are black. Furthermore, for every state S, denote by hs the probability that, 
starting at the state S, we eventually reach the full state (i.e. the state where all vertices are black). In the 
following we denote for simplicity the probability hs^ by /ifc, for every fc = 0, 1, . . . , n. Note that ho = and 
hi = hciique, since Sq is the empty state and Si is the state with only one black clique vertex. In order to 

compute an upper bound of /ii, we can set the value h " and the values hp. for every fc > 1 to their trivial 

upper bound 1. That is, we assume that the state S- 



, as well as all states F^, where fc > 1, are absorbing. 
After performing these relaxations (and eliminating self loops) , we obtain a Markov chain, whose state graph 



is illustrated in Figure 2(b) 
For any fc = 1, . . . , , " 

^ ' ' log n 



1 in this Markov chain, 



hk = akhk+i + Pkhk-i + 7fc 



(15) 



where 



Notice now by (fTB|) that 



oik 



Ik 



^ k(n—k) 



n \3 n J 



fc( 



1 -X- n-k ^ 
3 ' n / 



in-k _ 4 



(16) 



r ''"''+^^ +fc(H— ' 

n V, 3 n > 



> — > 1 



ak r{n — fc) 3r 



(17) 
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since r £ (l, |) by assumption. Furthermore, since i < | for sufficiently large n, it follows that for 



every k — 1,2,.. 



log n 



1, 



7fc 



1 7 

< 



Qffe n — k 6n 
Now, since afe + /3fe + 7fc = 1, (fT5|) implies by (fTT]) and (fT8| that 



^fc+l 



ak oik 



> iri^k 

6r 



1_ 



Thus, since /iq = and > hk-i for aU fc = 1, , 



i/c+i 



-) 

3rJ 

-) 



' log n ' 
k 



it fohows that for every fc. 



7 



A; 

) hi 



4_\k 

'-Sri 

6n' f - 1 

3r 



Consequently, since h " = 1 in the relaxed Markov chain, we have that 



(18) 



1-hi = 



{hk+i " hk) 

k=l 



> 



E 

fc=i 



/ii - — • ^ 

3r 



fe=0 



< 1 



6n 



y log' n 

3r fc=0 



and thus 



7 



1 



This completes the proof of the lemma, since > 1. ■ 
The next corollary follows by the proof of Lemma [TJ 

Corollary 1 Starting with one black clique vertex, the probability that we make at least one ring vertex black, 
or that we eventually reach ^^^y ^ black clique vertices, is at most J — + o (i). 

In the remainder of this section, we will also provide an upper bound on /i^mg , thus bounding the fixation 
probability /g„ of Gn (cf. Theorem HJ . Consider the Markov chain that is depicted in Figure [21 Our 
analysis will use the following auxiliary lemma which concerns the expected time to absorption of this Markov 
chain. 



P — p 
Figure 3: The Markov chain A4 



Lemma 2 Let p ^ q and p + q = 1. Then, as m tends to infinity, the expected number of steps needed for 
Ai to reach v„i, given that we start at vi, satisfies 

( ^™lnf+o(™) ^^^^^ 
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Proof. For i = 0, 1, . . . , m, let fii denote the expected number of steps needed to reach Vm- Clearly, /i^ — 
and /io = 1 + /ii- Furthermore, for i = l,...m — 1, it follows that 



^li = 1 +pm+i + qi-H-i 



(20) 



I.e. 



Consequently, we have that 



q, ^ 1 

P P 




(21) 



7+1 ~ Mi 



-Ml ^ 



Ml 



Ml 



m — 1 

E 

i=l 



1 



q-p 
1 



m 



S - 1 
p 



P 



which concludes the proof of the lemma for large m. 



(22) 



Denote in the following by A^i the Markov chain of the stochastic process defined in [T^ (see Section [TT] 
for an overview), when the underlying graph is the clique- wheel Gn, cf. Figure 2(a) The next definition will 
be useful for the discussion below. 

Definition 1 (Ring steps) A transition of the Markov chain Aii is called a ring step if it results in a 
change of the number of black vertices in the outer ring (i.e. ring vertices). 

We now present some domination statements that simplify the Markov chain A4i. More specifically, all 
these statements will increase the probability of reaching fixation when we start with one black ring vertex, 
such that we finally get an upper bound on hrmg- 

Di: Let w be a vertex on the outer ring, and let w' be its (unique) neighbor in the clique. We will forbid 
transitions of the Markov chain Mi, where a white colored v' tries to make white a black colored v. 

Z?2: Fixation is forced when either of the following happens: 
All The outer ring reaches logn black vertices. 

A2: The number of ring steps in order to reach logn black ring vertices is more than 0(log^ n). 
A^: The clique reaches n black vertices. 

^4: A black clique vertex makes black a white vertex on the outer ring. 

Let now be the modified Markov chain after these domination statements are imposed. Given these 
statements, we can prove the following. 

Proposition 1 If M2 ho,s not yet reached an absorbing state, then the outer ring consists of a single arc of 
consecutive black vertices, or has no black vertex at all. 

Proof. This follows by noticing that a second black arc can only be created either by a white colored vertex 
on the clique making white its black neighbor on the outer ring (this is impossible because of -Di), or by a 
black vertex on the clique making black its white neighbor on the outer ring (this will lead to absorption at 
Ai of D2). ■ 



The following definitions will be useful in the sequence. 
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Definition 2 (Off"spring) Given the history Ht of the Markov chain up to time t, the current state S QV 
at time t (i.e. the set of black colored vertices in the graph), and a vertex v € S , we will say that v is an 
offspring of u € V if and only if there exists a transition path in Ht that proves that v is black because of u. 

Notice in Definition [2] that u is not necessarily black at time t. 

Definition 3 (Birth in the clique) We will say that a vertex v' is born in the clique if and only if its 
(unique) neighbor v in the outer ring is black and makes a transition to the clique. 

Notice in Definition |3] that the color of v' is irrelevant before v' is born in the clique. We only need that 
the color of v is black. Furthermore, the above definition allows for a specific vertex to be born more than 
once (i.e. at different time steps). The proof of our main theorem can now be reduced to a collection of 
lemmas. Lemma [3] concerns the behavior of the ring. 

Lemma 3 Let Bi be the stochastic process describing the ring steps in Markov chain Ai2- Given that we do 
not have absorption at A4, then Bi is a birth- death process with forward bias equal to r. Furthermore, given 
that we start with a single black vertex on the ring, the following hold: 

(1) The probability that the number of black vertices in the outer ring reaches \ogn before absorption at 
A2,A3 or A4 is at most - — ^-j-y^F^- 

(2) The probability that more than log^ n ring steps are needed in order to reach logn black colored vertices 



in the ring, or to reach absorption in A2, A3, or A4 is at most O ^j^^tj^ ■ 



Proof. Recall that we do not allow transitions where the clique affects the number of black colored vertices 
in the outer ring (by the domination statements Di and A4). Then, it can be easily seen that the forward 
bias of the birth-death process Bi (i.e. the ratio of the forward probability over the backward probability) is 

2r 1 

-^-f = r, where W is the sum of the fitness of every vertex in the graph. Thus, part (1) of the lemma follows 

W 3 

by equation Q (for an overview of birth-death processes, see also pOll?!]). 

For part (2), let X denote the number of ring steps needed in order to reach logrt black colored vertices in 
the ring, or to reach absorption in A2, A3 or A4. Then X is stochastically dominated by the number of steps 
needed for Markov chain (cf. Figure [3]) to reach with m — logn and p = -p^. Hence, by Lemma [5] 
and Markov's inequality, we get that 

Pi{X > log^ n) <c-^ 



logn 

for some positive constant c = c(r). ■ 

The next lemma bounds the number of vertices that are born in the clique (see Definition [3]) . 

Lemma 4 Given that we start with a single black vertex on the ring, the probability that we have more than 
log^ n births in the clique is at most O „ ^ ■ 

Proof. For the proof, we will ignore for the moment what happens in the clique and how the clique affects 
the ring, since these steps are either forbidden (by Di) or lead to absorption (by A4). 

Let Y be the number of births in the clique (see Definition [3]) that we observe between two ring steps. 
Notice that at any time before absorption, there will be exactly 2 white colored vertices in the outer ring 
that can perform a ring step (see Definition [T]). Furthermore, if the number of black vertices in the ring is 
more than 2, then not all black vertices can affect the number of black vertices in the ring. We now restrict 
ourselves, to observe only ring- involved moves (forgetting about the clique), that is, transitions where only 
vertices of the ring that can cause a ring step or a birth in the clique are chosen. Given that M2 (i.e. the 
modified Markov chain) has not been absorbed, the probability that a ring step happens next is 

_ 2(1 + r) 1 

Pstep ^ r> I o 

2 + zr 3 

where z is the number of black colored vertices in the outer ring. Similarly, the probability that a birth in 
the clique happens next is 

zr 1 

Pbirth = n T7 

2 + zr S 
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Consequently, the random variable y + 1 is stochastically dominated by a geometric random variable with 
probability of success 

Pstep _ 2r + 2 ^1 



P = 



Pstep + Pbirth zr + 2r + 2 log n 
where in the last inequality we used the observation that at any time before absorption, the number of black 
vertices in the ring is at most log n because of Ai . But then, by Markov's inequality, we have that 

i 1 

Pr(F + 1 > log^ n + 1) < 5^ < ^ 

log n + 1 log n 

But by part (2) of Lemma [31 the probability that there are more than log^'n births in the clique before 
the Markov chain is absorbed is by Boole's inequality at most 

log^ nPr(r > log^ n) + o( <o(-^ 

\\ognJ \logn^ 

which concludes the proof of the lemma. ■ 

The following lemma states that it is highly unlikely that the clique will affect the outer ring, or that the 
number of black vertices in the clique will reach n. 

Lemma 5 Given that we start with a single black vertex on the ring, the probability of absorption at A3 
or A/^ is at most O ^1^^777^ • 

Proof. For the purposes of the proof, we assign to each birth in the clique a distinct label. Notice that, by 
Lemma m we will use at most log^ n labels with probability at least 1 — (^is^tt) • using more 

than log^ n labels (which happens with probability at most O (^j^;^^ by Lemma [4|), then we stop the process 

and assume that we have reached one of the absorbing states. Furthermore, whenever a black vertex v in 
the clique with label i replaces one of its neighbors with an offspring, then the label of v is inherited to its 
offspring. 

In order for M2 to reach absorption at A3, the clique must have n black vertices. Since each of these 
vertices has a label j G [log''n], there exists at least one label i such that at least y—r^ vertices have label i. 
Similarly, if reaches absorption at A4 and v is the corresponding affected ring vertex, then there exists a 
label i, such that v has label i. We will call a label i winner if there are at least ^ vertices in the clique 
that have label i, or the outer ring is affected by a clique vertex of label i. Clearly, if A^2 reaches absorption 
at ^3 of A4, there must be at least one winner. 

Recall that, by Corollary [1] the probability that a single black vertex in the clique either reaches j^^^-^ 

offsprings or affects the outer ring is at most ^ ^ J — + o(i). Consider now a particular label i. Then, 

if all the other black vertices of the graph that do not have label i (i.e. black ring vertices or black clique 
vertices with label j ^ i) had fitness 1, then the probability that i becomes a winner is by Corollary [1] at 
most — — \- + o (—). The fact that the other black vertices that do not have label i have fitness r can only 

6n(^-l) ^ri' 

reduce the probability that i becomes a winner. Therefore, considering all different labels i £ [log^ n] and 
using Boole's inequality, we conclude that the probability to reach absorption at ^3 or A4 is at most 




log' ni ' ^^ +0 i +0 -^^O 



log n j \ log 71 ^ 

where the term O {^j;^^ in the left side corresponds to the probability that we have more than log'' n labels. 
■ 

Finally, the following theorem concerns the probability of absorption of A^2- 

Theorem 3 For n large, given that we start with a single black vertex on the ring, the probability that M2 
is absorbed at Ai is at most (1 + o(l)) (l — i). Furthermore, the probability of absorption at A2, A3, or A4 

IS at most O (j^i^^) . 

Proof. The bounds on the absorption at Ai or A2 follow from Lemma [31 while the bounds on absorption at 
A3 or A4 follow from Lemma [SI ■ 

Recall now that A42 (the modified Markov chain) dominates A4i (the original Markov chain). Further- 
more, recall that the clique-wheel graph G„ has n clique vertices and n ring vertices, and thus the fixation 
probability of Gn is /g„ = ^'"'"'"''2^^'''"'' ■ Therefore, the next theorem is implied by Theorem [3] and Lemma [T] 
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Theorem 4 For the Markov chain Mi, and any r G (l, |), h^mg < (l ^ 7) • Therefore, as n 00, 



the fixation probability of the clique-wheel graph Gn in Figure 2(a) is by Lemma{l\ 



/G„<i(l-i)+o(l) (23) 

4 Analysis of the aggregation model 

In this section, we provide analytic results on the new evolutionary model of mutual influences. More 
specifically, in Section 14.11 we prove that this model admits a potential function for arbitrary undirected 
graphs and arbitrary initial fitness vector, which implies that the corresponding dynamic system converges 
to a stable state. Furthermore, in Section |4?2] we prove fast convergence of the dynamic system for the case 
of a complete graph, as well as we provide almost tight upper and lower bounds on the limit fitness, to which 
the system converges. 

4.1 Potential and convergence in general undirected graphs 

In the following theorem we prove convergence of the new model of mutual infiuences using a potential 
function. 

Theorem 5 Let G — {V^E) be a connected undirected graph. Let r(0) be an initial fitness vector of G, and 
let rmin and rmax be the smallest and the greatest initial fitness in r(0), respectively. Then, in the model of 
mutual influences, the fitness vector r{k) converges to a vector [rg, tq, . . . , tq]"^ as k 00, for some value 

Proof. Denote the vertices oi G by V ~ {ui,U2, . . ■ ,Un}. Let A: > 0. Then ^ implies that for any 
i = 1,2, ... ,n, the element r^,. (fc + 1) of the vector r{k + 1) is 

S fc ^ ' deg(uo) 

and thus 

r„,(fc + l) _ r,,(fc) 1 ^ r.,(fc)-r„,(fc) 

degK) degK) S(fc) ,^^.^^^^, "^-^ ^ deg(u.)deg(u,) ^' 

Therefore, by summing up the equations in (|24p for every i = l,2,...,nit follows that 

r„,(fc + l) ^ ru, (fc) ^ y- [ru^jk) - ruAk)f 

"^^^("0 ^^^deg(u,) S(fc)^^^^ deg(u,)deg(u,) 

> ru. (fc) 

Define now the potential function (/)(fc) = Tlu ev dc^lu'-) every iteration fc > of the process. Note by 
Observation [2] that S(fc) = X^u gy^"i(^) — "''max is a trivial upper bound for S(fc). Therefore, (I25p implies 
that 

. 1 V- (r„,(fc)-r„,(fc))^ 1 V- f ,,..2 
- A f \A t \ > ~3 ('•u,(fc) - r„,(fc)) 

"''max ^„ deg(Mi) deg(Uj ) 71'^rmax ^„ 

Furthermore, note that rmax ' X^m ev dcg(u ) — "''max is a trivial upper bound for (/)(fc). Therefore, since 
(j){k + 1) > 0(fc) for every fc > by (pS]) . it follows that ^(fc) converges to some value (^o as fc — ^ 00, where 
0(0) < 00 !i 'iT'max- Consider now an arbitrary e > and let e' = — . Then, since 0(fc) — > (po, there 

" nnax A; -1.00 
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exists ko € N, such that \(j){k + 1) — 4>{k)\ < e' for every k > fco- Therefore, (|26|) imphes that for every edge 
UiUj G of G and for every k > ko, 

K.(fc)-r„,(A:))2 < J2 i^u^k)-r.^^{k)f (27) 

< nVmax-|<^(A: + l)-0(A:)| < nV^ax • e' = 

Thus, for every e > 0, there exists fco G N, such that |r„^(fc) — rui{k)\ < e for every fc > fco and for every 
edge UiUj G £' of G. Therefore, since G is assumed to be connected, ah vahies ru{k), where u ^ V, converge 
to the same value tq as fc — >■ oo. Furthermore, since r„(fc) G [rmin , Tmax] by Observation [21 it foUows that 
''o G [?'min , ''max] as welL This completes the proof of the theorem. ■ 



4.2 Analysis of the complete graph 

The next theorem provides an almost tight analysis for the limit fitness value tq and the convergence time 
to this value, in the case of a complete graph (i.e. a homogeneous population). 

Theorem 6 Let G = (V,E) be the complete graph with n vertices and e > 0. Let a G [0, 1] be the portion 
of initially introduced mutants with relative fitness r > 1 in G, and let rg be the limit fitness of G. Then 
\ru{k) — ry{k)\ < e for every u,v €V, when 

fc > (n - 2) • ln(^— 



Furthermore, for the limit fitness rp , 

■ 1 + rv(r — 1 "1 4- 

1 + a(r - 1) 

and 



ro < 1 + a(r - 1) + ■ ^ ^ ' (28) 



> 1 + «(r - 1) + V(l + a{r - 1))^ + 2a(l - a)(r - 1^ ^29) 

> l + a(r-l) 

Proof. Since G is symmetric, we do not distinguish among the different placements S oi the an initially 
introduced mutants. Furthermore, at every iteration fc > 0, there exist by symmetry two different fitnesses 
ri(fc) and r2{k) for the vertices of S and of \ S', respectively. Thus, it suffices to compute only ri(fc) and 
r2(fc) for every fc > 0. Let A(fc) = ri(fc) - r2(fc). Then, A(0) = r - 1. It follows now by © and (0) that for 
every fc > 

(1 ~ a)nr2(k) ^ (1 — a)nr2(k) , ^ 

(1 — a)nr2{k) 



ri(fc) - A(fc) 



{n- l)S(fc) 



Similarly, 



where S(fc) — anri{k) + (1 — a)nr2{k). Subtracting now (|3T|) from ([30]), it follows that 

E(fc) 



A(fc + 1) = A(fc)-A(fc) 
= A(fc) 



(n- 1) • E(fc) 

n-2 



n — 1 

and thus, since A(0) = r — 1, it follows that for every fc > 



2\ k 



A{k) = {r-l)-[—j] (32) 
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Therefore, in particular, A{k) > for every fc > if and only if r > 1. Let now e > be arbitrary. Then 
|A(/c)| < £ if and only if 



n-2\'' e 

< r ^ 



n — 1 / r — 1 

r - 1 

> — (33) 



1 \'' r-1 



n-2 



1 + ■;^] — 7> e as n — >■ oo. Thus, for sufficiently large n, ([SS]) is satisfied when > or 
equivalently when 

k>{n~2)- ln(^^^) (34) 

Recall by Theorem [5] that ri(fc) — )> rg and r2(fc) — > vq for some value tq, as fc — )■ oo, and thus also 
ari(fc) + (1 — a)r2{k) — >• tq as fc — > oo. Furthermore, it follows by (|5n| and (PT|) that 

anik + 1) + (1 - a)r2(fc + 1) = «ri(fc) + (1 - a)r,ik) + "\ ,,,, • ^ (35) 

That is, ari(fc) + (1 — a)r2{k) is a non-decreasing function of fc, and thus ari(fc) + (1 — a)r2(fc) > ar + (1 — a). 
Therefore, for every fc > 0, 

ari(fc) + (1 - a)r2{k) < 1 + a(r - 1) + ^ "(^ ~ ") . J_ y A^k) (36) 



k=0 



The sum X^fc^o ^^(fc) can be computed by as 

i:A=W.(.-lf^.^-^.(,.-l)^(|-lI (37) 



Substituting now ([BT)) into ([BG)) . it follows that 



a(l — a) ^ 1 \2 " ^ 1 



ari(fc) + (1 - a)r2(fc) < 1 + «(r - 1) + — ^ • (r - 1)^- (38) 

1 + a[r — 1) 2n — 3 

Therefore, since — > i as n — >■ 00, and since ari(fc) + (1 — a)r2{k) — tq as fc — >■ 00, it follows by ([55)) 
that for sufficiently large n and fc, 

ro < 1 + a r - 1 -f ^ . \. ■ ^—^ 39 
1 + a(r — 1) 2 

Recall by ([55)) that Q;ri(fc) + (1 — Q:)r2(fc) is non-decreasing on fc, and thus ari(fc) + (1 — Q;)r2(fc) < tq. 
Therefore, it follows by (155)) and ([57)) that for every fc > 0, 

Q;(l — a) \2 ~ 1 



ari(fc) + (1 - a)r2(fc) > 1 + a(r - 1) H ^ ■ (r - I) 



ro 2n — 3 

Thus, since — > ^ as n — > 00 and Q;ri(fc) + (1 — Q;)r2(fc) — >■ tq as fc — ^ 00, it follows similarly to the above 
that for sufficiently large n and fc, 

ro > 1 + « r - 1 + ^ ^ • ^ „ ' 

ro 2 

and thus 

o , , XX Q;(l — a)(r — 1)^ 
r2-ro(l + a(r-l))-^ ^ ^>0 (40) 

Therefore, since ro > 0, it follows by solving the trinomial in ([^U)) that 



^ l + a{r-l) + V(l + a{r - 1))^ + 2a(l - a){F^ 
ro > ^ (41) 



The statement of the theorem follows now by ((34]), ([59)) . and (|4T 
The next corollary follows by Theorem [HI 
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Corollary 2 Let G = {V, E) be the complete graph with n vertices. Suppose that initially exactly one mutant 
with relative fitness r > 1 is placed in G and let tq he the limit fitness of G. Then 1 + - — - < t-q < 1 + ^ ^ ■ 



Proof. Since we have initially one mutant, it follows that a = ^. Then, substituting this value of a in ([29l) . 
we obtain the lower bound > 1 + For the upper bound of tq, it follows by substituting a in pS)) that 



ro < 1 



1 



r — 1 (r — 1 ) 



n 



71 ^ 71 ^ 



r — 1 , n — 1 r — 1, 

n r + (n — 1) 2 

r — 1 , r — 1 , 



2n 



5 Invasion control mechanisms 

As stated in the introduction of this paper, our new evolutionary model of mutual influences can be used to 
model control mechanisms over invading populations in networks. We demonstrate this by presenting two 
alternative scenarios in Sections 15.11 and 15.21 In both considered scenarios, we assume that an individuals 
of relative fitness r (the rest being of fitness 1) are introduced in the complete graph with n vertices. Then, 
as the process evolves, we periodically choose (arbitrarily) a small fraction j3 e [0, 1] of individuals in the 
current population and we reduce their current fitnesses to a value that is considered to correspond to the 
healthy state of the system (without loss of generality, this value in our setting is 1). In the remainder of this 
section, we call these modified individuals as "stabilizers" , as they help the population resist to the invasion 
of the mutants. 



5.1 Control of invasion in phases 

In the first scenario of controlling the invasion of advantageous mutants in networks, we insert stabilizers to 
the population in phases, as follows. In each phase k> 1, we let the process evolve until all fitnesses {r^ \ v G 
V} become e- relatively-close to their fixed point r^'^'' (i.e. until they e-approximate Tq'^-'). That is, until 
" — - < £ for every v . Note by Theorem[5]that, at every phase, the fitness values always e-approximate 

such a limit fitness r^\ After the end of each phase, we introduce fin stabilizers, where (3 G [0, 1]. That is, 
we replace fin vertices (arbitrarily chosen) by individuals of fitness 1, i.e. by resident individuals. Clearly, the 
more the number of phases, the closer the fixed point at the end of each phase will be to 1. In the following 
theorem we bound the number of phases needed until the system stabilizes, i.e. until the fitness of every 
vertex becomes sufficiently close to 1. 

Theorem 7 Let G = (V, E) he the complete graph with n vertices. Let a G [0, 1] be the portion of initially 
introduced mutants with relative fitness r > 1 in G and let j3 e [0, 1] he the portion of the stahilizers introduced 
at every phase. Let r^^ he the limit fitness after phase k and let £,S > 0, he such that ^ > ^/e and S > ^\/£. 
Finally, let each phase k run until the fitnesses e-approximate their fixed point r'^\ Then, after 

fc>l+ , (43) 
phases, the relative fitness of every vertex u is at most 1 + 5. 

Proof. Consider the first phase, where initially there exist an mutants with relative fitness r and (1 — a)n 
resident individuals with fitness 1 each. Then, since r > 1, it follows by for the fixed point rg^-* after the 
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first phase that 



(1) ^ 1 , ^ 1^ n , (1 - - 1) 
r-Q < l + a(r-l)-(lH 



2(l + a(r-l))' 



2 1 + a(r — 1) 

a(r~l) , 1, 
2 a 



I.e. 

r^'^<l + ^(r-l) (45) 

Suppose that we let each phase fc > 1 run until the fitnesses e-approximate their fixed point r^^K Note that, 
at the start of the process, (1 — a)n vertices have fitness 1 and an vertices have fitness r. Similarly, before 
the fcth phase starts, /3n vertices have fitness 1 and {l — /3)n vertices have fitness at most (1 + £)?'q'^~^''. Then, 
we obtain similarly to (|45p that the fixed point at iteration k is in the worst case 

Therefore 

(1 + < (1 + £) + (1 + e){l ^)((1 + e)r^t'^ - 1) 

and thus 

(1 + e)rf - 1 < e + (1 + - ^)((1 + e)rt'^ ~ 1) 
Let now A = (1 + £)(! — 'I). Then the last inequality becomes 

(1 + £)r(*=^ - 1 < e + A((l + e)rll'-''^ - 1) 

and by induction we have 

fe-2 

{l + e)/^^^! < e^A'- + A'=-i((l + £)r|,') -1) 

1=0 

1 - A'=-i 



1- A 

Therefore, (|45l) implies that 



+ A'^-n(l + £)rr-l) 



(1 + £)rr - 1 < ^^Y^ + + + 1)) (46) 

At the end of the fcth phase, the relative fitness of each vertex is at most (l + £)r'Q'^''. Now, in order to compute 
at least how many phases are needed to reach a relative fitness (1 + £)?'q'^'' < 1 + (5 for every vertex u € V, it 
suffices by (|46p to compute the smallest value of k, such that 

^^T^ + ^'-'(e + {l + s)^ir~l))<S (47) 

Recall now that ^/e < f < 5 by assumption. Therefore A = (1 + £)(! — ^) < (1 + £)(! — y/e), i.e. A < 1. 
Thus 1 — A*^~^ < 1 and it suffices from (jTf| to compute the smallest number k for which 

^ + A'^-i(£ + (1 + £)ii^(r -1))<S (48) 

Note now that 



1-A i_(i + e)(l-f) 

£ 

!(!+£)-£ 
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Thus, since f > s/e by assumption, it foUows that 



' < r." . = r (49) 



1-A Ve(l + e)-e 1 + e - Ve 

However 1 + e — ^/e > | for every e G (0, 1), and thus it foUows by (|49p that < \\fs. Therefore it suffices 
from (j48p to compute the smallest number k for which 

^Vi+A^-i(£ + (l+£)i^(r-l))<^. 



That is. 



e + il + e)^{r-l) 
or equivalently 



This completes the proof of the theorem. ■ 



5.2 Continuous control of invasion 

In this section we present another variation of controlling the invasion of advantageous mutants, using our 
new evolutionary model. In this variation, we do not proceed in phases; we rather introduce at every single 
iteration of the process /3n stabilizers, where /3 G [0, 1] is a small portion of the individuals of the population. 
For simplicity of the presentation, we assume that at every iteration the /3n stabilizers with relative fitness 1 
are the same. 

Theorem 8 Let G = {V, E) he the complete graph with n vertices. Let a G [0, 1] be the portion of initially 
introduced mutants with relative fitness r >\ in G and let j3 G [0, 1] he the portion of the stahilizers introduced 
at every iteration. Then, for every (5 > 0, after 

fc> ^(n-l).ln(^) (50) 

iterations, the relative fitness of every vertex u ^ V is at most \ -\- 5. 

Proof. Recall that we assumed for simplicity reasons that at every iteration the fin individuals with relative 
fitness 1 are the same. Note furthermore that at very iteration k we have by symmetry three different fitnesses 
on the vertices: (a) the an initial mutants with fitness ri(fc), (b) the jSn "stabilizers" with fitness 1, and (c) 
the rest (1 — a — /3)n individuals with fitness r2{k), where 1 < r2{k) < ri{k) by Observation [5] Note that 
''2(0) = 1. Let 7 = 1 — a — /3. Then, we obtain similarly to ([50]) and (|?T|) in the proof of Theorem [S] that for 
every A; > 

= ri{k) - (;^3Y)^(T"'^2(fc)(ri(fc) - r^ik)) + ^nik) - 1)) 

and 

,N ari(k)n ,,, (ari (fc) + /3)n , , /3n . 

'^^'^'^ - in - Dm • "^'^'^ + ' (n - i)m ^ - "^^^'^ + in - mk) ^''^ 

= r2ik) + l^(anri(fc)(ri(fc) - r2(fc)) - /3n(r2(fc) - 1)) 

where = n(ari(fc) + 7r2(fc) + /3). It follows now by ([SI]) and ([52]) that 

ri(fc + l) -r2(fc + l) = ri(fc)-r2(fc) 

(anri(fc) + 7nr2(fc))(ri(fc) - r2(fc)) + ^n(ri(fc) - r2(fc)) 



'"i(^) - r2{k) - 



(n- l)S(fc) 
S(fc)(ri(fc)-r2(fc)) 
in - l)I](fc) 



18 



and thus 

n-2 

ri{k + 1) - r2(fc + 1) = -(ri(fc) - r^ik)) 

n — 1 

Therefore, since r2(0) = 1 and »'i(0) = r > 1, it folfows that for every fc > 0, 



ri(fc)-r2(fc) = (r-l). ( — ^ ) (53) 



By substitution of into ((5T|) it foUows that 

ri(fc + l)=ri(fc)-^— ^^(7r2(A:)(r-l) (^-— + /3(ri (fc) - 1)) (54) 
Define now A(fc) = ri(fc) - 1. Then, it fohows by that 

< A(fc).(l ^ 



{n- l)S(fc)- 

^ fL n^A /f^^ 



Note now that > ^, and thus ([55]) imphes that 



A(fc + l)<A(fc). (1--^-^) (56) 

n-l- 



. /3 



Denote now for the purposes of the proof A = 1 — — n^x' ■ Then, it fohows by the system of 

inequahties in that for every fc > 

A(fc) < A(0) • A'= (57) 
= (r - 1) • A'^ 

In order to compute at least how many iterations are needed such that ri(fc) < 1 + (5, i.e. A(fc) < 5, it suffices 
by ([57|) to compute the smahest value of fc, such that 

(r - 1) • A*^ < (5 



I.e. 



1 , n — 1 , J. r — 1 , , 

r 

(1 + I )^- > '-^ 



However, (1 + xtttz^t^) ^ e as n — > oo. Thus, for sufficiently large n, ((55|) is satisfied when 



r- 1 
S 



or equivalently when 



r V — 1 

k > -^{n - 1) • ln(- 



This completes the proof of the theorem. ■ 

Observation 3 The hound in Theorem o/ the number of iterations needed to achieve everywhere a suffi- 
ciently small relative fitness is independent of the portion a G [0,1] of initially placed mutants in the graph. 
Instead, it depends only on the initial relative fitness r of the mutants and on the portion (3 € [0,1] of the 
vertices, to which we introduce the stabilizers. 
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6 Concluding remarks 



In this paper we investigated alternative models for evolutionary dynamics on graphs. In particular, we first 
considered the evolutionary model proposed in |16) , where vertices of the graph correspond to individuals of 
the population. We provided in this model generic upper and lower bounds of the fixation probability on 
a general graph G and we presented the first class of undirected graphs (called clique-wheels) that act as 
suppressors of selection. Specifically, we proved that the fixation probability of the clique-wheel graphs is 
at most one half of the fixation probability of the complete graph (i.e. the homogeneous population) as the 
number of vertices increases. An interesting open question in this model is whether there exist functions /i(r) 
and /2 (r) (independent of the size of the input graph) , such that the fixation probability of every undirected 
graph G lies between /i(r) and f2{'r)- Another line of future research is to investigate the behavior of the 
model of [16] in the case where there are more than two types of individuals {aggressive vs. non-aggressive) 
in the graph. 

As our main contribution, we introduced in this paper a new evolutionary model based on mutual in- 
fluences between individuals. In contrast to the model presented in |16| . in this new model all individuals 
interact simultaneously and the result is a compromise between aggressive and non-aggressive individuals. In 
other words, the behavior of the individuals in our new model and in the model of [16j can be interpreted as 
an "aggregation" vs. an "all-or-nothing" strategy, respectively. We prove that our new evolutionary model 
admits a potential function, which guarantees the convergence of the system for any graph topology and any 
initial fitnesses on the vertices of the underlying graph. Furthermore, we provide almost tight bounds on the 
limit fitness for the case of a complete graph, as well as a bound on the number of steps needed to approx- 
imate the stable state. Finally, our new model appears to be useful also in the abstract modeling of new 
control mechanisms over invading populations in networks. As an example, we demonstrated its usefulness 
by analyzing the behavior of two alternative control approaches. Many interesting open questions lie ahead 
in our new model. For instance, what is the speed of convergence and what is the limit fitness in arbitrary 
undirected graphs? What happens if many types of individuals simultaneously interact at every iteration? 

Acknowledgment. Paul G. Spirakis wishes to thank Josep Diaz, Leslie Ann Goldberg, and Maria Serna, 
for many inspiring discussions on the model of |16) . 
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